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广义 Vandermonde 方 程 组 的 有 效 快速 算法 * 
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(西北 工业 大 学 应 用 数学 系 ， 西 安 710072) 


摘 ”要 : 基于 求解 Vandermonde 方 程 组 的 Bjorck-Pereyra 算法， 本 文 给 出 了 求解 广义 Vandermonde 方 
程 组 的 有 效 快速 算法 ， 所 需 的 计算 量 为 O(n?)。 数 值 算 例 表明 ， 与 求解 Vandermonde 方 程 组 
的 Gohberg-Kailath-Koltracht 算法 和 Gauss 消 元 法 相 比 ， 本 文 的 算法 具有 更 高 的 计算 精度 。 
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1 引言 


定义 aaz ,an 是 nn 个 实数 或 复数 ， 称 如 下 的 nw 阶 方 阵 
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AJ X Vandermonde 矩阵 。 

许多 实际 的 问题 可 以 转化 为 广义 Vandermonde 和 矩阵 的 相关 求解 问题 ， 如 要 构造 次 数 不 超 
过 nn 的 缺 项 多 项 式 g(x) = co + cix + coz? +--+ ent"? + cuz" EnS A a, a2, ,an 处 
满足 插值 条 件 g(akr) = f&(k = 1,2,-,n) 3X — Inl UE (JR E] X. Vandermonde 7j f£ 
组 Vie = f， 其 中 c = (co ,cn 20) = (fifo o Sn) o WRR- RDE, 
如 Gauss 消 元 法 ， 求 解 广义 Vandermonde 方 程 组 ， 所 需 运 算 量 为 O("3)。 近 年 来 ， 针 对 一 
些 特殊 矩阵 出 现 了 一 批 求解 相应 计算 问题 的 快速 算法 ， 其 运算 量 达到 O(n2) RERS, 
对 于 由 Vandermonde 和 矩阵 4 = (oj ' ?ij-1 推 广 得 到 的 所 谓 Vandermonde 型 矩阵 (本 文 的 广 
义 Vandermonde 和 矩阵 是 其 特例 )， 文 献 2 利 用 这 类 和 矩 阵 的 位 移 结构 给 出 了 求解 相应 方程 组 
的 快速 算法 : Gohberg-Kailath-Koltracht 算法 ， 但 由 于 Vandermonde (型 ) 矩阵 一 般 是 严重 病 
态 的 ， 用 所 给 快速 算法 和 Gauss 消 元 法 求解 相应 方程 组 的 阶 数 很 低 。 文 献 [1] 利用 Newton 插 
值 多 项 式 构造 了 求解 Vandermonde 方 程 组 的 有 效 快速 算法 : Bjorck-Pereyra 算 法 。 本 文通 
EK n I X Vandermonde EREA RN A n 十 1 阶 的 Vandermonde 和 矩阵， 并 利用 分 块 处 理 技巧 
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和 求解 Vandermonde 方 程 组 的 Bjorck-Pereyra 算法， 构造 出 了 求解 广义 Vandermonde 方 程 组 
的 快速 算法 ， 所 需 的 运算 量 为 O(n?)。 数 值 算 例 表明 ， 与 Gauss 消 元 法 和 Gohberg-Kailath- 
Koltracht 算法 相 比 ， 本 文 的 算法 具有 更 高 的 计算 精度 。 


2 ”快速 算法 推导 
考虑 求解 广义 Vandermonde 方 程 组 
V'z-f, Vy-g (2) 
HP V ÆR (1) K X Vandermonde &BPE,. Tfj 


£ = (Z1 £n)", y = (yi, Yn) fee g (9,94). 


HE n + 18t Vandermonde 4E RE 


t7 P n+l 
V = à; hagan (3) 


EP ai, a2, Qn, Qm 十 1 两 两 互 异 ， Banı # 0. WP- (el ,en_1) en+ly En)» 其 中 ei 是 
第 i 个 分 量 为 1， 其 余 分 量 为 0 的 n 十 1 维 列 向 量 。 则 有 


其 中 
Ga (iia a) A Vo ryan lo) Tod. 
令 
"20 * Des. 
0T 1 
结合 式 (4) 可 得 


上 式 两 边 求 道 ， 并 利用 Sherman-Morrison 公式 和 分 块 矩阵 求 道 公式 得 


1 A a m y- -Vla 
人 i 
[s] 0 -IEV- Vota) 
0 


pi 
s= (51,77 vnd t= (tis ,tnt1)’, © = (Wi, ;Un41)^, z= (nv uad 
分 别 是 Vandermonde 方程 组 


m-(1) Vi=enn, Vé-g Vi-h (6) 
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的 解 向 量 ， 其 中 


D T rd 二 T 
g= (gj, , gn—1, 0, ga ) > h= (1,a 1,7: an 0,08 a) $ 


romuna (1) wa tyme (! 


9 外 auraa 
0 0 


s= (si. ， Sm 一 1， 3Sm 十 1 s t= (te ,tn—1, tn 二 1 ys 


w = (ww Wn), z= (222,77 28) - 


则 由 式 (7) 和 (8) 可 得 广义 Vandermonde 方程 组 (2) 的 解 为 


T 


as 
g-—sd4———i y—-wvwct 


Un-4-1 
i-ant E (9) 


工 一 Zn+1 


式 (6) 的 Vandermonde 方 程 组 可 直接 由 Bjorck-Pereyra 算 法 求解 ， 综 上 可 得 如 下 的 快速 算法 。 
算法 1 KN X Vandermonde HFEA VIr = f. 


i—1 。 = nn . 
ai = ap; (i= 1,2,- n- 1), Qn = ana 
十 


: 1 1 
c = fis pt =0 (j m 12, n), ct， = 0, Pb =1. 


(k) — (k) (k) 0 
ders Cj 二 Cj-1 poe Pi Pj- 


了 7 了 Deep 2 
Qj Qj—k 


j-kctlkc2,-,n4l, 
Qj 一 Qj—k 


十 1 (n1 +1) _ +1 
d = cn+1 k E kz pet ). 
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对 k= 二 n,.… ,2,1， 有 


一 OkSy 11 , j=k+1, 7 


= 


k k+1 1 1 . 

QD QD sysh, 4D, j=1,2,. n1 
n—1 n—1 

o= Y ois; 十 ansn+l，， HS 》 aiti 十 antn+l， Zn = 9841 Ola a /(1 — p), 
i—1 i—l 


ri = si +ot;/(1— pu), i—1,2,--,n-— 1. 


该 算法 需 (5n?/2) + (9n/2) + 工 次 乘除 运算 和 3n2 + 7n 一 2 次 加 减 运算 。 
算法 2 求解 广义 Vandermonde 方 程 组 Vy = g. 


1 1 j— 。 
d? = gi, gP mal, j—712,-.,n—-1, 


1 1 n 
d(? 一 gt) — 0, d, = n; a. = Qni: 


对 k= 1,2,- ,D5 有 


= a(? m: adi?) uen E i is arq, j= k+1,k +2, ,n+1, 


1 j 1 j ; 
wt Lr E dr di j71,2,-,ncl 


对 k= 二 n,.… ,2,1， 有 


ufk+1) E (b+ ， ute+1) = P i—1,2,--. k, 


T zc wt /(a; — 83k), - ind zm x7 f (aj " jk); j= k+ lk +2, n+ 1, 


w = D 


l 49 2 D, 


, 


i-21,2,-.-,k— 1, 


vf o cui, HP ouem aem, joker 

k k+1 k k+1 1 1 : 
up sut E Mores a wj = uf ? aja ), J21L2,.- ,n- 1, 
Yi = Wi + Wagizi/l(l—zn41) i=1,2, ,n. 


该 算法 需 (5n2/2) + (5n/2) -- 1 次 乘除 运算 和 3n2? + 5n 次 加 减 运算 。 


3 ”数值 算 例 


我 们 针对 一 些 广 义 Vandermonde 和 矩阵 V， 在 微机 (CPU 奔 IW2.6G， 内 存 256MB) 上 编程 
计算 分别 用 本 文 所 给 的 快速 算法 、Gauss 消 元 法 和 Gohberg-Kailath-Koltracht 算法 ( 简 记 
为 GKK FA) RH X Vandermonde JHA Vz = f 和 Vy = g， 部 分 算 例如 下 。 

例 1 RH X Vandermonde FFEA VIz = f, HP 


ai =—1 +2(i—1)/n, fi =Tn-1lai) +27? (a? a?) i=1,2,-,n, 
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这 里 Ti(z) 是 m — 1r Chebyshev 多项式。 为 方便 起 见 ， 取 nn 为 奇数 。 该 方程 组 的 精确 解 
为 Chebyshev E MA 了 T,_1(x) 的 展开 式 系数 ， 即 


ape E DEEE a (RR, a= (CARB 
DG 


Wan = 1， 误 差 用 e = ||z* 一 zllz 度 量 ， 计 算 结果 见 表 1。 


表 1: 广义 Vandermonde 方 程 组 VTz = f 算 例 结果 比较 


15 19 23 27 
2.74e-08 | 3.24e-06 6.01e-04 | 9.30e-01 


3.22e-06 | 3.46e-03 6.88e+00 | 
4.12e-01 | 3.16e+07 


n T 11 
快速 算法 1 1.77e-13 | 8.13e-11 
高 斯 消 元 法 | 3.86613 | 1.16e-09 
GKK 算 法 5.84e-11 | 2.62e-06 


| 


例 2 求解 广义 Vandermonde 方 程 组 Vy = g, H 


= Ji? i=1,2,.……,n—1, n = 


这 里 cf 二 n(n 一].…(n 一 7 上 DAL。 该 方程 组 的 精确 解 为 
yj = (7157 0(1-3/2" ^, j=1,2, n. 


H anı = EE 误差 用 
ly; — vil 
e= max 一 一 一 
iSj*" — [yjl 


度量 ， 计 算 结果 见 表 2. 


表 2: J À Vandermonde 2 4248 Vy = g 算 例 结果 比较 


16 | 20 
1.76e-10 | 7.91e-06 | 3.65e-01 
2.27e-05 | 1.03e+00 

1.00e--00 


n 4 8 
快速 算法 2 | 9.47e-15 | 8.27e-13 
高 斯 消 元 法 | 1.38e-14 | 1.25e-09 
GKK 算 法 | 5.92e-16 | 1.54e-05 


数值 算 例 表明 ， 与 Gauss 消 元 法 和 Gohberg-Kailath-Koltracht 算法 相 比 ， 本 文 的 算法 具有 
更 高 的 计算 精度 。 由 于 广义 Vandermonde 和 矩阵 是 严重 病态 的 ， 所 以 以 上 算 例 中 和 矩阵 的 阶 数 不 
高 。 本 文 研究 的 广义 Vandermonde 方 程 组 的 快速 算法 很 容易 推广 到 和 矩阵 为 一 般 跳 行 的 情形 。 
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An Efficient and Fast Algorithm for Solving Generalized 
Vandermonde Systems 


ZHAO Liang-dong, XU Zhong, LU Quan 


(Department of Applied Mathematics, Northwestern Polytechnical University, Xi'an 710072) 


Abstract: In this paper, an efficient and fast algorithm for solving generalized Vandermonde systems 
is obtained on basis of the Bjorck-Pereyra algorithm for solving Vandermonde systems. The algorithm 
costs O(n?) arithmetic operations. Numerical results show that the algorithm is higher in precision 
than the Gauss elimination and Gohberg-Kailath-Koltracht algorithm when solving Vandermonde-type 


systems. 
Keywords: Vandermonde matrix; generalized Vandermonde matrix; linear systems; fast algorithm 
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